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ABSTRACT 

The eclipsing and double-lined spectroscopic binary system V453 Cygni consists of two early 
B-type stars, one of which is nearing the terminal age main sequence and one which is roughly 
halfway through its main sequence lifetime. Accurate measurements of the masses and radii of 
the two stars are available, which makes a detailed abundance analysis both more interesting 
and more precise than for isolated stars. We have reconstructed the spectra of the individual 
components of V453 Cyg from the observed composite spectra using the technique of spectral 
disentangling. From these disentangled spectra we have obtained improved effective tempera- 
ture measurements of 27 900 ± 400 K and 26 200 ± 500 K, for the primary and secondary stars 
respectively, by fitting non-LTE theoretical line profiles to the hydrogen Balmer lines. Armed 
with these high-precision effective temperatures and the accurately known surface gravities 
of the stars we have obtained the abundances of helium and metallic elements. A detailed 
abundance analysis of the primary star shows a normal (solar) helium abundance if the micro- 
turbulence velocity derived from metallic lines is used. The elemental abundances show no 
indication that CNO-processed material is present in the photosphere of this high-mass ter- 
minal age main sequence star The elemental abundances of the secondary star were derived 
by differential study against a template spectrum of a star with similar characteristics. Both 
the primary and secondary components display elemental abundances which are in the ranges 
observed in the Galactic OB stars. 
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1 INTRODUCTION 

High-mass stars are major producers of ultraviolet light and 
byproducts of thermonuclear fusion, so are important objects in the 
study of star formation, and galactic chemical and kinetic evolu- 
tion (e.g. Samland 1998). The large luminosities of high-mass stars 
means that they are important distance indicators in both a galac- 
tic and cosmological context (e.g. Guinan et al. 1998; Ribas et al. 
2005; Southworth et al. 2005, 2007). 

In the last decade new theoretical models have been con- 
structed which incorporate a number of physical processes impor- 
tant to the structure and evolution of high-mass stars, including 
convective core overshooting, semiconvection, rotational mixing 
and mass loss. Overshooting and rotation can both have a large 
effect on the predicted lifetimes and luminosities of high-mass 
stars (c.f. review by Maeder & Meynet 2000). However, empiri- 
cal constraints on these processes remain hard to come by, despite 
a steady improvement in observational techniques and capabilities 
(see Hilditch 2004). 

Detached eclipsing binaries (dEBs) are vital objects for ob- 
taining observational constraints on the structure and evolution of 
high-mass stars, since they are the primary source of directly mea- 
sured stellar properties (Andersen 1991). Unfortunately, accurate 



(2% or better) physical properties are available for only eight high- 
mass and most of these objects have no observational con- 
straints on their chemical composition. Chemical abundances are 
difficult to determine from the spectra of high-mass dEBs for sev- 
eral reasons. Firstly, they tend to display only a small number of 
spectral lines. Secondly, the often high rotational velocities of these 
objects means their spectral lines are wide and shallow so high 
signal-to-noise ratio (S/N) spectra are needed to obtain useful re- 
sults. Thirdly, in dEBs the spectral lines from one component inter- 
fere with the line profiles of the other component ('line blending'). 

In a seminal work, Simon & Sturm (1994) introduced the tech- 
nique of spectral disentangling (SPD), by which individual spec- 
tra of the component stars of a double-lined spectroscopic binary 
system can be deduced from a set of composite spectra observed 
over a range of orbital phases. SPD can be used to measure spec- 
troscopic orbits which are not affected by the line blending which 
afflicts other methods, such as cross-correlation (see Southworth & 
Clausen 2007). The resulting disentangled spectra also have a much 
higher S/N than the original observations, making them very useful 



^ An up-to-date compilation of the properties of well-studied dEBs is avail- 
able at http : / / www .astro . keele ■ ac ■ uk/~ jkt| 



2 K. Pavlovski and J. Southworth 



for chemical abundance analyses. As a bonus, the strong degener- 
acy between effective temperature (Tcs) and surface gravity (log g) 
is not a problem for dEBs because surface gravities can be mea- 
sured to high accuracy (0.01 dex or better) for both components. A 
detailed study of these possibilities is given by Pavlovski & Hens- 
berge (2005; hereafter PH05). 

1.1 The eclipsing binary system V453 Cyg 

V453 Cyg belongs to rather sparse group of high-mass close bi- 
nary systems for which accurate physical properties have been mea- 
sured for both components. In the most recent study of V453 Cyg, 
Southworth, Maxted & Smalley (2004b; hereafter SMS04) anal- 
ysed new spectroscopic observations and previously published 
light curves, obtaining masses of Ma ~ 14.36 ± 0.20 Mq and 
Mb = 11.11 ±0.13 Mq, and radii of i?A = 8.55 ± 0.06 Rq and 
Rb = 5.49 ± 0.06 Kq. The high precision of the radius measure- 
ments was helped by the geometry of the V453 Cyg system, which 
exhibits deep and total eclipses. SMS04 derived effective temper- 
atures of Tcft A = 26 500 ± 800 K and Teff b = 25 300 ± 600 K 
from the equivalent widths of several He I and Hell lines, using 
an analysis based on local thermodynamic equilibrium (LTE) syn- 
thetic spectra but with a correction for non-LTE effects. 

V453 Cyg is known to display the phenomenon of apsidal mo- 
tion, which can be used as a probe of the internal structure of the 
components of an eccentric binary system (e.g. Claret & Gimenez 
1993). SMS04 found an apsidalperiod of t/ = 66.4± 1.8 yr, which 
gives the structural parameteo logfc2 — —2.254 ± 0.024. The 
masses, radii, Tcff s and log k2 of V453 Cyg match theoretical pre- 
dictions for an age of 10.0 ± 0.2 Myr and an approximately so- 
lar metal abundance. The primary component (the hotter and more 
massive star) is reaching the terminal age of its main sequence life- 
time (TAMS) whilst the secondary star is roughly halfway through 
its MS life. The significant difference between the properties of the 
two stars means that V453 Cyg is a very useful system for testing 
evolutionary models. 

In this work we analyse the spectroscopy of V453 Cyg ob- 
tained by SMS04 and by Simon & Sturm (1994) in order to ob- 
tain new TcffS and accurate chemical abundances for the stars. 
V453 Cyg was targeted by SMS04 as part of a program to study 
dEBs which are members of open clusters (see Southworth, Maxted 
& Smalley 2004a; Southworth et al. 2004c). It is a probable mem- 
ber of NGC 687 1 , so our results can be compared to published anal- 
yses of single stars which are also cluster members. 



2 SPECTROSCOPIC DATA 
2.1 Calar Alto spectra 

A total of 29 spectra of V453 Cyg were obtained by Simon & 
Sturm (1994) during two observing runs in 1990 and 1992 with 
the 2.2 m telescope at the German-Spanish Astronomical Centre 
at Calar Alto, Spain (Table[T](. Four different spectral intervals 
were observed with the coude spectrograph, each covering 250 A. 
The intervals are centred on the hydrogen lines US (A4050), H7 
(A4430), H/3 (A4850), and on a 4533-4730 A region containing 
metallic lines. Hell A4686 and He I A4712. The reciprocal disper- 
sion is 0.13 A px~^ for the 1990 spectra and 0.20 A px~^ for the 

^ Note that the log ^2 value in the abstract of SMS04 is incorrect: the cor- 
rect value is given in SMS04 (their Section 7.2) and here. 



Table 1. Observing log for the Calar Alto spectra of V453Cyg. Ac is 
the central wavelength and the orbital cycles are calculated from HJD 
2453 940.0998. 



Set 


ID 


HJD 


S/N 


Ac (A) 


Cycle 


A 


21x 


48813.39567 


150 


4010 


2435.4036 


A 


43x 


48814.63754 


150 


4010 


2435.7229 


A 


51x 


48816.38687 


160 


4010 


2436.1726 


A 


61x 


48817.38987 


170 


4010 


2436.4304 


A 


71x 


48818.37820 


190 


4010 


2436.6845 


A 


81x 


48819.38231 


190 


4010 


2436.9427 


A 


32 


49201.45130 


140 


4010 


2535.1653 


A 


41 


49202.35881 


170 


4010 


2535.3986 


A 


51 


49203.46754 


190 


4010 


2535.6837 


B 


22x 


48813.47704 


180 


4290 


2435.4245 


B 


42x 


48815.46905 


170 


4290 


2435.9366 


B 


62x 


48816.52872 


220 


4290 


2436.2091 


B 


72x 


48817.64888 


200 


4290 


2436.4970 


B 


85x 


48818.60773 


190 


4290 


2436.7435 


B 


82x 


48819.47156 


190 


4290 


2436.9656 


C 


23x 


48812.57408 


180 


4520 


2435.1924 


C 


32x 


48814.48183 


180 


4520 


2435.6828 


C 


41x 


48815.40237 


140 


4520 


2435.9195 


C 


52x 


48816.42184 


190 


4520 


2436.1816 


C 


63x 


48816.63435 


210 


4520 


2436.2362 


C 


83x 


48818.59129 


200 


4520 


2436.7393 


D 


llx 


48812.48692 


200 


4750 


2435.1700 


D 


11 


49199.49495 


200 


4750 


2534.6624 


D 


12 


49199.60984 


190 


4750 


2534.6919 


D 


13 


49199.65734 


150 


4750 


2534.7041 


D 


21 


49200.38908 


170 


4750 


2534.8923 


D 


31 


49201.41599 


180 


4750 


2535.1562 


D 


33 


49201.66650 


150 


4750 


2535.2207 


D 


42 


49202.38781 


170 


4750 


2535.4061 



1992 data. The S/N ratios are 140-220; further details are given by 
Simon etal. (1994). 



2.2 La Palma spectra 

A set of 43 spectra of V453 Cyg were secured in 2001 October 
with the 2.5 m Isaac Newton Telescope and Intermediate Disper- 
sion Spectrograph at La Palma, Spain (Table|2j- These are split 
between two spectral intervals of length 260 A, reciprocal disper- 
sion 0.11 Apx^^ and resolution 0.20 A. Most spectra are centred 
on 4360 A and the remainder on 4580 A: both datasets cover the 
He I A447 1 and Mg II A448 1 lines, which are good lines for radial 
velocity (RV) determination over a wide range of spectral types 
(Andersen 1975). A few additional spectra were obtained with cen- 
tral wavelengths of 4340 A and 4840 A; they are not used here. The 
S/N of the La Palma spectra is 100-450 per pixel, with most close 
to 150. We have rejected a few spectra with a S/N below 100. 



3 METHOD 

Our aim is to perform a detailed test of the predictions of stel- 
lar evolution models of high-mass stars using strict empirical con- 
straints. This aim can only be reached by studying dEBs, as these 
are the only objects where we can measure masses and radii of 
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Table 2. Observing log of La Palma spectra of V453 Cyg. Columns are the 
same as for Table[T] 



Set 


ID 


HJD 


S/N 


Ac (A) 


Cycle 


E 


323747 


52562.50016 


160 


4360 


3399.2270 


E 


323867 


52563.31850 


110 


4360 


3399.4374 


E 


324029 


52564.31988 


160 


4360 


3399.6948 


E 


324255 


52565.33945 


180 


4360 


3399.9569 


E 


324455 


52566.29608 


140 


4360 


3400.2029 


E 


324462 


52566.31116 


140 


4360 


3400.2068 


E 


324505 


52566.43565 


180 


4360 


3400.2388 


E 


324581 


52568.34709 


160 


4360 


3400.7302 


E 


324757 


52569.31004 


100 


4360 


3400.9777 


E 


324774 


52569.34031 


130 


4360 


3400.9855 


F 


323085 


52559.33775 


210 


4580 


3398.4140 


F 


323086 


52559.34009 


210 


4580 


3398.4146 


F 


323285 


52560.31872 


150 


4580 


3398.6662 


F 


323286 


52560.32104 


150 


4580 


3398.6668 


F 


323287 


52560.32335 


150 


4580 


3398.6674 


F 


323288 


52560.32567 


150 


4580 


3398.6680 


F 


323289 


52560.32799 


150 


4580 


3398.6686 


F 


323317 


52560.40911 


130 


4580 


3398.6894 


F 


323318 


52560.41142 


130 


4580 


3398.6900 


F 


323319 


52560.41375 


130 


4580 


3398.6906 


F 


323320 


52560.41606 


130 


4580 


3398.6912 


F 


323321 


52560.41837 


130 


4580 


3398.6918 


F 


323479 


52561.30395 


140 


4580 


3398.9195 


F 


323480 


52561.30627 


150 


4580 


3398.9201 


F 


323481 


52561.30858 


150 


4580 


3398.9207 


F 


323482 


52561.31091 


140 


4580 


3398.9213 


F 


323483 


52561.31322 


150 


4580 


3398.9219 


F 


323744 


52562.49486 


160 


4580 


3399.2257 


F 


323864 


52563.31216 


110 


4580 


3399.4358 


F 


324026 


52564.31111 


170 


4580 


3399.6926 


F 


324247 


52565.32035 


180 


4580 


3399.9520 


F 


324584 


52568.35585 


170 


4580 


3400.7324 


F 


325129 


52570.35446 


170 


4580 


3401.2462 


F 


325348 


52571.30720 


180 


4580 


3401.4911 



high-mass stars to accuracies of 1%. In this context the SPD tech- 
nique has an important role both in mass and abundance determi- 
nation. RV and mass measurements of OB-type dEBs are difficult 
because their spectra generally contain only a few lines, which have 
strong rotational broadening and so suffer line blending. 

Using SPD and a set of spectra taken at a range of orbital 
phases, it is possible to determine the individual component spec- 
tra. One can also simultaneously solve for the RVs for each compo- 
nent and spectrum, or directly for the orbital parameters of the sys- 
tem. Contrary to the cross-correlation method (Simkin 1974; Tonry 
& Davis 1979) no template spectrum is needed since the individual 
component spectra are themselves used as templates, so a signifi- 
cant source of bias is bypassed (c.f. Hensberge & Pavlovski 2007). 

Several studies have put extensive effort into reconciling RVs 
measured by several techniques to orbital parameters derived using 
SPD. Orbital parameters from RVs determined by multi-component 
Gaussian fitting to complex line profiles have been found to be infe- 
rior to those derived by SPD (Simon & Sturm 1994b; Simon, Sturm 
& Fiedler 1994; Hensberge et al. 2000). SPD also performs better 
than standard cross-correlation (Holmgren et al. 1997). 

Southworth & Clausen (2007) critiqued these methods and 
the two-dimensional cross-correlation technique TODCOR (Zucker 
& Mazeh 1994) in an analysis of the dEB DWCarinae. The main 
problem in RV determination was found to be line blending. SPD 



Table 3. Parameters of the spectroscopic orbits derived in this work and by 
SMS04. SMS04 fixed the orbital period (P = 3.889825 d), time of pe- 
riastron passage (Tperi = 2 439 340.6765) and eccentricity (e) to values 
found in their apsidal motion analysis. In the present work we have opti- 
mised e, the longitude of periastron (a;) and the amplitudes of RV variation 
(Ka and K^). 



Parameter 


SMS04 


Subset C 


Subset F 


e 


0.022 (fixed) 


0.029 ±0.016 


0.024 ± 0.003 




140.1 (fixed) 


153.4 ± 13.7 


166.7 ±8.6 


Ka (kms-i) 


173.7 ±0.8 


174.1 ±0.9 


172.5 ±0.2 


i^B(kms-i) 


224.6 ±2.0 


216.2 ± 1.7 


221.5 ±0.5 


q 


0.773 ±0.008 


0.805 ± 0.008 


0.779 ±0.003 



gave the most reliable RVs, and fitting a double-Gaussian also gave 
acceptable results. Cross-correlation techniques could not properly 
account for line blending and gave results which were system- 
atically biased towards smaller stellar masses. These conclusions 
were obtained in a case when the two component stars had very 
similar properties and individual spectral lines were studied, so 
some of the benefits of cross-correlation could not be accessed. 



4 SPECTROSCOPIC ORBITS THROUGH SPECTRAL 
DISENTANGLING 

Two spectroscopic datasets are available for analysis which were 
obtained about 10 years apart (Section^. As V453 Cyg exhibits 
fast apsidal motion (see Section fTTH these cannot be analysed to- 
gether. The datasets are split over a number of spectral windows: 
subset F (La Palma) is the most extensive and includes 24 spec- 
tra centred on 4580 A. This subset was used for determination of 
the orbital parameters. The phase distribution of the spectra is rea- 
sonable but suffers slightly from the orbital period of V453 Cyg 
(3.89 d) being close to an integer number of days: phases = 0.0- 
0.20 could not be covered. However, an observation was deliber- 
ately obtained during secondary eclipse, where only the primary 
star is visible. This spectrum was actually used as the template in 
the TODCOR analysis performed by SMS04 as it is a good match 
to both stars. Disentangled spectra formally contain no information 
about the continuum level for each star (PH05), but in the case of 
V453 Cyg this secondary-eclipse spectrum can be used to reliably 
determine how much light each star contributes. 

We have performed a spectral disentangling analysis in veloc- 
ity space (see Hadrava 1995) using the FDBINAR'vQ code (Ilijic et 
al. 2004). In Fourier SPD individual pixel flux uncertainties cannot 
be used but individual spectra can be assigned a weight based on 
their quality. We weighted the input spectra according to their S/N 
and adopted as our initial orbit the final solution given by SMS04. 
We fixed the orbital period at 3.889825 d (SMS04). 

SPD was performed over the spectral region 4540-4690 A. The 
strong lines of He II A4471 and Mg II A4481 were excluded in the 
determination of the orbital parameters, as in our experience better 
results are obtained for larger numbers of weak lines than for a few 
strong lines. The resulting parameters are given in Table|3] The er- 
rors are determined by the jackknife method (e.g., Yang & Robin- 
son 1986), which seems to be an appropriate approach for error 
calculations in SPD. Table[3]also contains the final orbit of SMS04, 
which is in good agreement with the results found here. Note that, 

^ |http : / /sail ■ zpf ■ f er ■ hr /f dbinary/l 
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Figure 1. Example region of the observed composite spectra (red lines) from La Palma subset F. The two individual disentangled spectra are shown at the 
bottom of the plot (blue lines) with their correct continuum levels. The disentangled spectra have then been adjusted with the appropriate Doppler shifts and 
relative intensities to reproduce the observed spectra, and are overlaid on them using blue lines. All spectra have been normalised to the continuum, and the 
observed spectra have each been shifted by +0.1 from their neighbour for clarity. The spectrum taken during secondary eclipse is the sixth from the bottom 
among the observed spectra. 



in their analysis, SMS04 used TODCOR to measure RVs so were 
only able to use spectra where the two stars had very different RVs. 
They thus fixed the eccentricity, e, and periastron longitude, uo, to 
the values found from their apsidal motion analysis. In our solution 
we were able to use all the available 4580 A spectra (subset F), in- 
cluding that obtained during secondary eclipse (phase 0.50). The 
orbital parameters are therefore more reliable, particularly e and uo. 
An example of SPD of the observed composite spectra and resulting 
disentangled spectra are shown in Fig.[T] 

For the Calar Alto spectra, subset C was chosen to be the mas- 
ter dataset for determining the orbital parameters. Since the indi- 
vidual Calar Alto datasets contain fewer spectra than the La Palma 
ones, we fixed the eccentricity to the value obtained from subset F 
(e = 0.023). Other parameters were adjusted in SPD. Other subsets 
(A, B, and D) were disentangled keeping the orbital parameters 



fixed to those obtained for subset C. The La Palma spectra were 
disentangled using the orbital parameters obtained for subset F. 

A feature of SPD is that the disentangled spectra are effectively 
co-added from the original observations and so have a higher S/N. 
The total S/N in a disentangled spectrum of component n is 



Sdis ~ /riSobsViVobs 

where /„ is the fractional light contribution for component n, and 
Sobs and A^obs are the S/N and number of input spectra. Both sets 
of the La Palma spectra cover the 4455^505 A region (which in- 
cludes He I A4471 and Mgll A4481), so the total number of ob- 
servations of this spectral window is 45. We have disentangled this 
region separately, yielding individual component spectra with S/N 
of about 700 for the primary and 250 for the secondary (/a = 0.75 
and /b = 0.25). For comparison, disentangling the eight spectra 
in subset C (Calar Alto) gives individual spectra with S/N of about 
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380 and 130. The enhanced S/N in the disentangled spectra of the 
components relative to the observed spectra can be seen in Fig.[T] 
The two disentangled spectra have different S/N because the sec- 
ondary star is fainter than the primary. 



5 SPECTRAL ANALYSIS OF BOTH COMPONENTS 

The primary goal of this study is to measure the photospheric ele- 
mental abundances of the components of V453 Cyg, which requires 
knowledge of their Tcff S and loggs. Whilst disentangled spectra 
can be analysed in the same ways as single-star spectra, a potential 
drawback is that the continuum level for each star is not specified 
by SPD, short of ensuring individual spectral lines do not dip below 
zero flux. For V453 Cyg this is not a problem because spectra have 
been obtained during eclipses (spectrum 324777 in subset F and 
spectrum llx in subsetD). Also, the solution of the Cohen (1974) 
light curves of V453 Cyg (SMS04) allows us to assign a light factor 
to the two stars for each observed spectrum. We have renormalised 
the disentangled spectra using the method and formulae presented 
by PH05. The mid-eclipse spectra are themselves suitable for the 
following analysis, but we have used the renormalised disentangled 
spectra due to their much higher S/N. 

Abundance analyses of single stars can suffer from a strong 
correlation between the T^s and log g measurements from stellar 
spectra. In the case of dEBs this can be avoided because the log g 
of each star can be measured accurately. The values for V453 Cyg 
are log 3a = 3.731±0.012 and log = 4.005±0.015 (SMS04). 

5.1 Effective temperature determination 

SMS04 used the equivalent widths of helium lines and LTE the- 
oretical spectra (with a correction for non-LTE effects) to derive 
Toft A = 26 500 ± 800 K and Tofta = 25 300 ± 600 K, but noted 
that the helium abundance might be enhanced and that this could 
affect the measurements. As the helium abundance is an important 
part of the current analysis, we could not use helium lines to mea- 
sure TeffS. 

A standard way of deriving T^s for a sequence of B stars is 
from the Si Il/Si Ill/Si IV ionisation equilibrium (c.f. Becker & But- 
ler 1990). Some spectral lines of these three ions are visible in the 
primary's disentangled spectra, but its relatively high projected ro- 
tational velocity (llOkms"^) means the equivalent width of Sill 
A 4128 is difficult to measure to sufficient accuracy. Also, some 
important lines (e.g. Si II A5379) are not covered by the available 
spectra. In the secondary's spectrum. Si IV A 4212 is barely visible. 

We have instead used the Balmer lines, which are good indi- 
cators of Teff in hot stars if log g is already known from elsewhere. 
Our disentangled spectra cover H/9, H7 and US. In OB stars the 
Balmer lines are affected by numerous oxygen, nitrogen and silicon 
lines, so we carefully preselected 'clean' parts suitable for match- 
ing to theoretical line profiles. Grids of theoretical spectra were cal- 
culated for Tcff = 24 000-29 000 K and for the known log gs of the 
components of V453 Cyg. For these calculations we used ATLAS9 
model atmospheres (Kurucz 1974) with a microturbulence veloc- 
ity of nturb ~ 2kms^^ and solar metal abundances. The codes 
DETAIL (Giddings 1981; Butler 1994) and SURFACE (Butler 1984) 
were used for the non-LTE calculations - DETAIL determines the 
atomic level populations by jointly solving the radiative transfer 
and statistical equilibrium equations, and SURFACE computes the 
emergent fluxes and line profiles. Model atoms were used for H 
(Herrero 1987; Husfeld et al. 1989), He (Husfeld et al. 1989), C II 



(Eber & Butler 1988), Nil (Becker & Butler 1989), Oil (Becker 
& Butler 1988), Mgll (Przybilla et al. 2001), AlIII (Dufton et al. 
1986), and Si II, Si III and Si IV (Becker & Butler 1990). 

FORTRAN programs have been developed which fit theoretical 
line profile to hydrogen lines by minimisation in Teg. More 
weight was given to the H7 profile since it was derived using a 
dataset which contains a spectrum taken during mid-eclipse, and is 
covered by both the Calar Alto and the La Palma data sets. Parts of 
the profiles which are affected by metallic lines were excluded from 
the fit. The best fits were found for T^sa = 27 800 ± 400 K and 
TcftB = 26 200 ± 500 K, and are plotted in Fig.[2]for the H7 and 
H/9 lines. These TeffS are larger by 1. 5a and l.la than the values 
found by SMS04, which can be ascribed to the importance of non- 
LTE effects at these temperatures. Note that this does not affect the 
light or RV curve analyses presented by SMS04. 

While not used in deriving the above TeffS, the lines of He II 
were examined for the reliability of the results. Our data cover three 
He II lines, at at A4199, A4541 and A4686. A broad and diffuse line 
at A4199 is barely visible in the primary star's spectrum. The lines 
at A4541 and A4685 are consistent with derived primary T^b but are 
better fitted with T^a ~ 27 500 K (0.75cr lower than found above). 
The match of the Hell line at A4686 in the secondary's spectrum 
with theoretical profiles for T^b = 26 200 K is excellent. 

The projected rotational velocities (v sin i) of the components 
were derived from the widths of the clean spectral features due to 
CII A4267, Sim A4553, Sim A4568, Oil A4591, Oil A4596, 
and O II A4662, avoiding He I and Mg II lines (Hensberge et al. 
2000), using a set of the theoretical spectra calculated for dif- 
ferent vsini values. We find «Asini = 109 ± 3kms"^ and 
vb sin j = 98 ± 5kms~^. These figures are in excellent agree- 
ment with those found by SMS04. The pseudo-synchronous ro- 
tational velocities are llOkms"^ and 71kms~^; the primary is 
therefore pseudo-synchronised but the secondary is rotating more 
quickly than expected. 

In their abundance analysis of disentangled component spec- 
tra of the eclipsing and double-lined binary V578 Mon, a mem- 
ber of the Rosette nebula cluster NGC 2244, PH05 used a sharp- 
lined spectrum of the cluster member star #201 as a template spec- 
trum. Vrancken et al. (1997) derived T^g = 26 500 ± 1000 K and 
log (7 — 4.3 ± 0.2 for this star, so it should be a good match to 
V453CygB. A rotationally broadened spectrum of #201 closely 
resembles V453 Cyg B, which lends additional support to the TcftS 
derived in this work. 

As an independent check we have used Stromgren uvbyp cal- 
ibrations to obtain the mean Toff of the system. The b — y, mi 
and ci indices given by Crawford & Barnes (1974) and Reimann 
(1989) are in good agreement with each other, indicating that they 
are not affected by the photometric variability of V453 Cyg. The 
UVBYBETA code was used to apply the calibration of Moon & 
Dworetsky (1985) to the Stromgren indices. We found a redden- 
ing of Eb-y — 0.33 mag, which is in excellent agreement with the 
value of Eb-v = 0.46 ± 0.03 mag found by Massey, lohnson & 
DeGioia-Eastwood (1995) for NGC 6871, and a T^g of 27 720 K. 
This Tcft value is a flux-weighted mean for the two stars because 
we do not have individual Stromgren indices for the two compo- 
nents, and is in good agreement with the values we find from the 
Balmer line profiles. 
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Figure 2. Comparison between the disentangled spectra (points) and best-fitting theoretical spectra (solid lines) for H(5 (left panel), H7 (middle panel) and H/3 
(right panel). Spectra for the secondary star have been offset by —0.2 for clarity. 



5.1.1 Implications of the new effective temperature 
measurements 

Armed with the T^s values determined above, we have calculated 
the distance to V453 Cyg using theoretical bolometric corrections 
published by Bessell, Castelli & Plez (1998) and by Girardi et al. 
(2002). We adopted the spectroscopic orbits of the two stars from 
Section|4] and the photometric parameters from SMS04. The two 
sources of bolometric corrections are in almost perfect agreement, 
and lead to a \/-band distance of dv ~ 1720 ± 80 pc and a if -band 
value of (Ik = 1680 ± 30 pc. The uncertainty in the reddening 
towards V453 Cyg dominates the uncertainty in dv but is much 
less important for dK- This distance is not in good agreement with 
the value of 2140 ± 70 pc found for NGC6871 by Massey et al. 
(1995), but is consistent with other values in the literature. 

We have compared the properties of V453 Cyg with several 
sets of theoretical models, but cannot find a good match for any of 
the chemical compositions for which models are available. Inter- 
estingly, the best (i.e. least bad) match is obtained with the Claret 
(1995) models for a 'solar' metal abundance (Z — 0.02), a high 
fractional helium abundance of F = 0.38, and an age of 7.4 Myr. 
In this case the Tcfts and the radius of the primary star are correctly 
predicted, but the measured radius of the secondary star is a few 
sigma smaller than the model value. A new photometric analysis of 
V453 Cyg would help to investigate this discrepancy further. 



5.2 Helium abundance 

The spectra of early-B stars are rich in strong hydrogen and He I 
lines. For the hottest B stars. Hell lines become visible. The in- 
clusion of non-LTE effects is important when determining the he- 
lium abundance of such stars, but its importance varies for different 
lines. The DETAIL and SURFACE codes use a model atom calcu- 
lated including transitions up to n = 4 (Przybilla & Butler 2001), 
so the A4471 and A4922 helium lines are calculated precisely but 
the diffuse A4026 and A4388 lines are calculated less precisely. 
A4713 is also reliable as it comes from a transition to level n — 4. 
As discussed in detail by Lyubimkov et al. (2004) regarding line 
broadening theory, the A447 1 and A4922 lines are to be preferred 
for helium abundance determinations. 

There is growing evidence that, in the calculation of helium 
line profiles, «turb should be carefully investigated (Gies & Lam- 
bert 1992; McErlean, Lennon & Dufton 1999; Lyubimkov, Ros- 



Table 4. Helium abundances derived for the primary component of 
V453 Cyg with the microturbulent velocity as free and fixed parameter 



Line 


e(He) 


Vturb free 


e(He) 


vturb fixed 


4026 


0.096±0.006 


5±3 


0.081 ±0.008 


15 (fix) 


4388 


0.113±0.008 


7±2 


0.075±0.010 


15 (fix) 


4471 


0.152±0.006 


6±3 


0.103±0.006 


15 (fix) 


4713 


0.118±0.009 


7±2 






4921 


0.134±0.009 


8±2 


0.086±0.012 


15 (fix) 


mean 


0.123±0.017 


7±2 


0.086±0.012 


15 (fix) 



topchin & Lambert 2004). Therefore, an extended grid of theo- 
retical spectra were calculated for the Tcss and loggs found for 
the components of V453 Cyg and for helium abundances from 
ehc = 0.09-0.21 in steps of 0.03, and t;turb ~ 0-20 km s~^ in 
steps of 5 km s~^. 

First, we have used the approach of Lyubimkov et al. (2004) 
to simultaneously determine wturb and eho- A4471 and A4922 were 
used for determining ehc as they have only a slight dependence on 
''turb- The A4713 line has a strong dependence on Wturb so was used 
to constrain this parameter once eno was determined. As before, 
minimisation was used. The results are given in Table|4](columns 2 
and 3). We find an enhanced helium abundance of ehc = 0.123 ± 

0. 017 for the primary star. It is evident that a large scatter between 
different He I lines is present. The two largest values of the helium 
overabundance were obtained for the He I A4471 and A4921 lines, 

1. e. the lines Lyubimkov et al. (2004) favoured in determination of 
the helium abundance in their sample. Conversely, we find a normal 
abundance of ehc = 0.089 ± 0.003 for the secondary star. These 
results are discussed further in Section[6] 

Secondly, we have determined the helium abundances using 
the value of uturb obtained from metallic lines, in particular oxy- 
gen. In hot OB stars the O II lines are relatively strong and nu- 
merous, so are well suited to this purpose. Uturb was derived by 
requiring the abundances from individual O II lines to be inde- 
pendent of the line strengths. The equivalent widths were mea- 
sured for 22 O II lines which are not severely blended. We find 
fturb = 15 ± 2kms~^. 

It is well documented by previous studies that the fturb tends 
to be larger for lower surface gravities (e.g. Kilian 1992; Vrancken 
et al. 2000; Daflon et al. 2004 and references therein). In the sample 
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Figure 3. The spectrum of the secondary component (sohd hne) of V578 Mon derived by SPD (Hensberge et al. 2000) overplotted on the disentangled spectrum 
of the secondary component of V453 Cyg (filled symbols, this work). The stars have almost the same T^ff , and their log gs differ by only 0. 15 dex. It is clear 
that the stars have very similar abundance patterns. 



of nine slowly-rotating (3 Cephei stars recently analyzed by Morel 
et al. (2008), three stars have Wturb > 10kms~^, whilst the full 
range of Vturb is 1-14 km s"'^. All of these objects are in similar 
evolutionary state as V453CygA (logg ~ 3.75) so our derived 
fturb should be reliable. 

Theoretical He I line profiles were calculated for ene = 0.07- 
0.15 and uturb ~ 15kms~^. The Tofi and logg were taken 
from Section [5T1 The best fit by minimisation is a mean he- 
lium abundance which is considerably different from results ob- 
tained with ?;turb as a free parameter. The mean helium abun- 
dance in the primary's photosphere from the four He I lines is 
e(iife) = 0.086 ± 0.012, almost equal to the canonical solar he- 
lium abundance [eQ(He) = 0.089]. The importance of this result 
is discussed in Section|6] 

5.3 Metals 

5.3.1 Primary star 

Determination of the Uturb of the primary star also yielded an oxy- 
gen abundance of log e(0) = 8.42 ± 0.087. This estimate is ob- 
tained from 22 O II lines, and the quoted error is the rms around the 
mean. The abundances of the other elements were estimated from 
fitting synthetic spectra to the disentangled spectrum of the primary 
star. A grid of synthetic spectra were calculated for the parameters 
of the star and different elemental abundances. The uturb was fixed 
to 15 km s~^, the value found from fitting the star's oxygen lines. 

About ten lines were available for the estimates of the abun- 
dances for nitrogen and silicon. The majority of N lines are blended 
by O II lines, but show a good interagreement for an abundance 
loge(N) = 7.80 ± 0.08 (11 lines). By contrast, our measurement 
of the Si abundance (log e(Si) = 7.27±0.14, 9 lines) suffers from a 
larger scatter in results for individual lines. The strong Si III triplet 
at A4550 to A4580 give larger abundances than the weak Si III lines 
on the far blue wing of H/3. This could be due to an incorrect ?;turb 



or to contamination by the H/3 line. Interestingly, the analysis of (3 
Cephei stars by Morel et al. (2008), also shows the largest disper- 
sion in abundances for Si. This issue is also discussed by Trundle 
et al. (2008) who used silicon lines for deriving fturb- 

Only four carbon lines are present in our spectra: the strong 
CII A4267, the weaker C III A4186, and two C II lines at A4650- 
A4680 which are heavily blended with O II lines. Our estimated C 
abundance (Table|5j suffers from the small number of useful lines 
and also from the quality of the available model atom (Trundle et 
al. 2008; Nieva & Przybilla 2008). 

The only visible magnesium line in the spectrum of the pri- 
mary star is Mg II A4481, which is blended with Al III A4479. Alu- 
minium is also detectable through other two lines (AlIII A4512 
and A4532). The abundances of Mg and Al were derived by calcu- 
lating theoretical spectra fitting blended line profiles in the 4477- 
4484 A region. The best match was found for loge(Mg) = 7.58 
and log e(Al) = 6.10. The other two Al III lines give abundances 
consistent with this. 

Beside intrinsic (random) errors, several systematic errors af- 
fect the determination of abundances, such as uncertainties in the 
atmospheric parameters (Tcff, \ogg and utmb) and in fundamen- 
tal atomic data {gf values). Morel et al. (2006) and Trundle et al. 
(2007) have discussed these problems, and we have adopted their 
approach in estimating the overall uncertainties in our abundance 
measurements. We repeated our abundance analysis using synthetic 
spectra with atmospheric parameters modified by AT = lOOOK, 
and Auturb ~ 3km s~^. The uncertainty in the log gs of V453 Cyg 
is only 0.01 dex, so is negligible here. The overall uncertainties in 
the abundances for the primary star (Table[5]l include the line-to- 
line scatter and the effects of changes in Tcff and uturb- The error 
budget was calculated only for O, N and Si, for which sufficient 
lines are present in the primary star's spectrum. For the other ele- 
ments, a conservative upper limit on the total error of 0.20 dex is 
assigned. 
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Table 5. Abundances derived for the components of V453 Cyg 



Abundance 


Primary 


Secondary 


OB stars 


H/He 


0.094±0.015(4) 


0.098±0.011 


0.10 


loge(C) 


8.15±0.20 (2) 


8.30±0.12 


8.2 


loge(N) 


7.80±0.10(11) 


7.60±0.10 


7.6 


loge(O) 


8.58±0.12(22) 


8.55±0.05 


8.5 


loge(Mg) 


7.58±0.20(1) 


7.50±0.10 


7.4 


loge(Si) 


7.27±0.16(9) 


7.25±0.12 


7.2 


loge(Al) 


6.12±0.20(3) 




6.1 



5.3.2 Secondary star 

We have so far concentrated on the primary star, as its mass and 
evolutionary state are astrophysically very interesting. We now 
measure elemental abundances for the secondary star by differen- 
tial analysis. In the study of the dEB V578 Mon, PH05 devised a 
complete methodology for the determination of elemental abun- 
dances from disentangled spectra of high-mass stars. Special at- 
tention was given to the minimisation of systematic errors arising 
from continuum placement. The abundance determinations were 
made relative to rotationally broadened spectra of a slowly-rotating 
star which was otherwise similar to the components of the binary. 
A detailed examination of the spectra and comparison to template 
spectra made it possible to extract the blended line profiles. Here 
we use the same methodology and wavelength ranges as in PH05. 

V453CygB (Tcff = 26 200 ± 500 K) closely resembles 
V578 Mon B (T^s = 26 400 ± 400 K). Their surface gravities dif- 
fer by only 0. 15dex and their rotational velocities are similar, so 
a direct comparison of their spectra is possible. Fig.[3]is a plot of 
the A4635-A4657 spectral region for the two stars, which contains 
a complex blend of O II, N II, N III, C II, C III and Si IV lines. 

There is striking similarity of both spectra throughout the ob- 
served spectral windows. Some slight differences can be found, and 
it is almost certain that these are intrinsic to the stars rather than an 
artifact of SPD. For comparison synthetic spectra have been cal- 
culated for elemental abundances of an 'average' main-sequence 
B star; mean values of the abundances are log e(C) = 8.27, 
loge(N) = 7.62, loge(O) = 8.57, loge(Mg) = 7.48 and 
log e(Si) = 7.25 (Daflon, Cunha & Becker 2004). The spectra of 
both stars closely follow the characteristics of the synthetic spec- 
trum. The elemental abundances of V578MonB resemble 'aver- 
age' B-dwarf abundances in the quoted errors (PH05, their table 3). 

As V453 Cyg is a probable member of the open cluster 
NGC 6871, which is part of the Cyg OB 3 association, we can com- 
pare our results with the work of Daflon et al. (2001). These authors 
determined elemental abundances of OB stars in five OB associa- 
tions. Spectra of eights objects from CygOB3 were obtained, but 
final analyses were performed for only three sharp lined targetfl 

The results of Daflon et al. (2001) for the stars in CygOB3 
are: loge(C) = 8.15 ±0.23, loge(N) = 7.68 ± 0.19, log £(0) = 
8.54 ± 0.14, log e(Mg) = 7.63 ± 0.24, log e(Si) = 7.22 ± 0.11. 
Our quantities for both V453 Cyg A and V453 Cyg B are in excel- 
lent agreement with these figures. 



* V453Cyg was included in the study by Daflon et al. (2001), who found 
T(,fl = 29 900 K and logg = 4.53. The composite nature of the spectra 
was not accounted for, so these quantities are not reliable. 



6 DISCUSSION 

Due to the serious problems of line blending and high rotational 
velocities, abundances have been measured for only a few high- 
mass stars in double-lined binary systems. Prior to the invention 
of SPD (Simon & Sturm 1994) some helium abundances in high- 
mass binaries were determined from equivalent widths measured 
using Gaussian fitting (Lyubimkov 1998; Pavlovski 2004 and refer- 
ences therein). The use of SPD has allowed helium abundances to be 
measured for four close binaries: DHCep (Sturm & Simon 1994), 
Y Cyg (Simon et al. 1994), V578 Mon (Hensberge et al. 2000), and 
DWCar (Southworth & Clausen 2007). Preliminary results have 
also been announced for V380Cyg (Pavlovski et al. 2005), and 
V478 Cyg and CW Cep (Pavlovski & Tamajo 2007). Finally, Frey- 
hammer et al. (2001) have derived helium abundances for the com- 
ponents of CPD — 59°2628 by fitting theoretical spectra. 

Lyubimkov et al. (2004) found an enrichment of helium during 
the MS evolution of their sample of stars, although many theoret- 
ical evolutionary models (e.g. Schaller et al. 1992) predict no sur- 
face enrichment of helium or other CNO-cycle elements. However, 
more recent theoretical studies (Maeder & Meynet 2000; Heger & 
Langer 2000) have found that rotationally induced mixing can af- 
fect the chemical composition of the surface layer of high-mass 
stars even whilst they are on the MS. Due to CNO processing in the 
stellar core. He and N should be enriched and C and O depleted. 
The strength of these effects depends on initial rotational velocity, 
so observational constraints can be used to both calibrate theoreti- 
cal models and check their predictions (e.g. Venn et al. 2000). 

The observations of 100 Galactic B stars by Lyubimkov et al. 
(2000, 2004) have confirmed that helium overabundance is corre- 
lated with stellar age expressed as a fraction of its MS lifetime. 
Large spectroscopic surveys of OB stars in our Galaxy and the 
Large and Small Magellanic Clouds (Hunter et al. 2008) have fur- 
ther confirmed this effect at a range of metallicities. Unfortunately, 
investigations using single stars suffer from the low accuracy with 
which their masses, radii and ages can be measured. Further work 
must therefore concentrate on binary stars, where all these proper- 
ties can be measured to within a few percent (e.g. SMS04). 

The subsample of high-mass stars in the sample of OB stars 
studied by Lyubimkov et al. (2004) show helium enrichment during 
their MS lifetimes, which has been found previously for close bi- 
naries (see Lyubimkov 1998 and references therein). In two recent 
studies in which large samples of B stars were analysed, progres- 
sively increasing helium enrichment from the ZAMS to the TAMS 
was found (Lyubimkov et al. 2004, Huang & Gies 2006). Both stud- 
ies saw this effect over their whole mass ranges, with large enrich- 
ments for high-mass stars. Moreover, Huang & Gies noted that the 
effect is large among the faster rotators (see their fig. 12). Lyu- 
bimkov et al. derived utmb both from helium lines and from O II, 
while Huang & Gies (working with only a limited spectral range 
around H7) accepted utmb ~ 2kms~^ for all their sample af- 
ter considering the discussion of uturb in Lyubimkov et al. (2004). 
Therefore, the agreement achieved by these studies may depend on 
their use of a common uturb- 

So far, in the sample of OB binaries studied by SPD and cited 
above, no MS star has been found to have a helium enrichment. To 
this list we add now both components of V453 Cyg. The ages of the 
stars on the list lie between and 0.8 of their MS lifetimes, with 
V453 Cyg A contributing the upper limit. Possible explanations for 
the non-detection of any changes in elemental abundances are mod- 
erate rotational velocity and tidal effects between the components. 
In this context it is interesting that detailed abundance analyses of 
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about a dozen (3 Cephei stars (Morel et al. 2006) has given an av- 
erage helium abundance of e(He) = 0.087 ± 0.012, with a range 
from 0.075 to 0.100. These objects are slowly-rotating stars in ad- 
vanced phases of their MS lifetimes and do not show enhanced he- 
lium abundances. The TcS and log g of V453 Cyg A fall within the 
range of their sample, and this star also does not show enhanced 
helium. These two results are in excellent agreement and corrob- 
orate previous finding of Herrero et al. (1992), who found a very 
clear dependence of helium enrichment on rotational velocity, but 
not on fractional lifetime on the MS, for 0-type stars. 

In Section lS^ we determined the He abundance in 
V453CygA by two approaches: (i) simultaneous line-profile 
fitting allowing ?;turb to be a free parameter; (ii) fixing ?;turb to the 
value obtained from O II lines. When using a fixed ?;turb we found 
a normal helium abundance (e(He) — 0.086 ± 0.012), but when 
Wturb was a free parameter we found a helium overabundance 
(e(He) = 0.123 ± 0.017) and also erratic metal abundances. We 
therefore reject the helium enrichment scenario and caution that it 
is crucial to carefully consider the value of uturb in similar studies. 



7 SUMMARY 

Binary and multiple stellar systems are our primary source of fun- 
damental measurements of the basic physical properties of stars. 
Modern observational data and analytical tools allow us to deter- 
mine stellar masses and radii with accuracies of 1 % or better, which 
in turn yields surface gravity measurements accurate to a few hun- 
dredths of a dex. When combined with advanced spectroscopic and 
abundance analysis techniques it is now possible to perform de- 
tailed tests of the rotation and magnetic effects included in modern 
theoretical evolutionary models. 

In this work we have studied the eclipsing double-lined spec- 
troscopic binary V453 Cyg, a system which contains two stars with 
accurately measured masses and radii. V453 Cyg A is approaching 
the TAMS and V453 Cyg B has completed about half of its MS 
phase. These characteristics make V453 Cyg ideal for testing the 
MS chemical evolution of high-mass stars, which is affected by 
convection, rotational mixing, and magnetic fields. 

Through SPD of several sets of phase-resolved high-resolution 
spectroscopy we have obtained the individual spectra of the two 
stars. By fitting the hydrogen line profiles with non-LTE synthetic 
spectra, we have measured the effective temperatures for the two 
stars to be T^sj^ = 27 800 ± 400 K and Teff b = 26 200 ± 500 K. 
The high precision of these measurements was achieved because 
the surface gravities of the two stars are known to within 0.01 dex 
(SMS04). The Tcff s are larger than have been found previously, and 
are supported by the Stromgren photometric indices of the system. 

The new Tcff values imply a distance of 1680 ± 30 pc to 
V453 Cyg (using theoretical bolometric corrections), which is in 
good agreement with most literature values but still shorter than 
the distance of 2 140 ± 70 pc found for NGC 687 1 by Massey et al. 
(1995). The predictions of theoretical stellar evolutionary models 
are unable to match the properties of V453 Cyg. 

Using the disentangled spectra we have measured the helium 
and metal abundances in the photospheres of the two stars by fitting 
non-UTE synthetic spectra to the observed line profiles. The metal 
abundances for the secondary star were estimated by differential 
analysis to a template star of similar Tefj and log g. We find that 
the helium in the photosphere of the primary has a solar abundance 
if we adopt a microturbulent velocity obtained from metallic lines. 
This result does not contradict evolutionary calculations (Meynet 



& Maeder 2000; Heger & Langer 2000) because a slight helium 
enrichment is predicted only for highly rotating stars, and/or after 
TAMS. Also, as is expected for the stars which were formed from 
the same interstellar material, no significant difference in the chem- 
ical composition between the two components is apparent. Their 
photospheric abundances resemble those of the typical Galactic OB 
stars (c.f. Herrero 2003 and references therein). 

We are planning to extend our analysis to larger numbers of 
spectral lines, which will require new spectroscopic observations 
with a wider wavelength coverage. This will make it possible to 
refine the present abundance analysis, and to address more thor- 
oughly the question of which Uturb to adopt for different ions. 
Echelle spectroscopic observations will be well suited to this task. 
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